Electron-Electron Interactions in Artificial Graphene 
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Recent advances in the creation and modulation of graphene- like systems are introducing a science 
of "designer Dirac materials" . In its original definition, artificial graphene is a man-made nanostruc- 
ture that consists of identical potential wells (quantum dots) arranged in a adjustable honeycomb 
lattice in the two-dimensional electron gas. As our ability to control the quality of artificial graphene 
samples improves, so grows the need for an accurate theory of its electronic properties, including 
the effects of electron-electron interactions. Here we determine those effects on the band structure 
and on the emergence of Dirac points. 

PACS numbers: 73.21.Fg, 73. 21. La, 73.22.Pr 



Graphene, a single two-dimensional (2D) crystal of car- 
bon atoms, has become the most attractive carbon-based 
material and one of the hottest topics in condensed mat- 
ter and material physics [1]. Many of the superlatives 
attributed to graphene are due to the fact that, at low 
energies, electrons (and holes) behave as massless chiral 
Dirac fermions as a result of a linear dispersion relation 
near two inequivalent corners of the Brillouin zone, i.e., 
the Dirac points. 

Not surprisingly, material combinations and arrange- 
ments with properties similar to real graphene have 
been sought [2, [3| , and experimental progress in produc- 
ing artificial structures has been impressive. In 2009, 
the use of a hexagonal nanopatterned 2D electron gas 
(2DEG) - similar to previously considered triangular an- 
tidot arrays 0, [H[ - was suggested [6] , and soon there- 
after Gibertini et al. Q presented the graphene-like band 
structure of GaAs quantum dots (QDs) arranged in a 
honeycomb lattice - a system they called as "artificial 
graphene" (AG). Their tight-binding calculations were 
supplemented by an experimental demonstration of a 
nanopatterned modulation-doped sample [7|, |8| . Very re- 
cently, the same authors and their collaborators showed 
that AG subjected to a strong magnetic field exhibits col- 
lective modes according to the Mott-Hubbard model [9j. 
In a parallel development, a structure equivalent to AG 
has been fabricated by STM-controlled deposition of CO 
molecules on the 2DEG on Cu (111) surface [3]. This 
leads to extremely controlled samples having even less 
defects than natural graphene. The trend limited to elec- 
tronic systems: as another way to control Dirac fermions 
Tarruell et al. [2] have created a tunable honeycomb op- 
tical lattice in an ultracold quantum gas. 

One of the reasons for pursuing the study of AG is 
that this system offers the opportunity to experimen- 
tally study regimes that are difficult or impossible to 
achieve in natural graphene, (including high magnetic 
fluxes, tunable lattice constants, precise manipulation of 



defects, edges, and strain) and allow the experimental ob- 
servation of several predictions made for massless Dirac 
fermions j!ol - [l2| . In Ref. [l3[ experimental criteria for the 
realization of graphene-like physics in 2DEG have been 
described. Recently, the system has been proposed as a 
candidate for the observation of a quantized anomalous 
Hall insulator [14j . 

These rapid experimental advances call for state-of- 
the-art numerical tools to calculate electronic properties 
of AG and, more generally, of artificial lattices. In view 
of the accuracy with which the sample can be prepared, 
it is particularly important to be able to reliably predict 
the conditions under which isolated Dirac points will ap- 
pear in AG. By "isolated Dirac points" we mean a set of 
points in momentum space where a conical conduction 
band makes contact with a conical valence band - the 
contact occurring at an energy at which no other state 
exists. As shown in Ref. Q such points occur only if the 
depth of the potential within the QDs of the AG struc- 
ture exceeds a certain minimum value. However, this 
minimum value cannot be reliably predicted from a the- 
ory that neglects the electron-electron (e-e) interactions. 
It is not just a matter of replacing the bare potential by 
an effective one that includes interaction effects such as 
screening, exchange and correlation. The key point is 
that this effective potential must be consistent with an 
electronic density distribution that places the Fermi level 
at the Dirac point. For example, the effective potential 
at N = 1 electrons per dot is quite different from the 
effective potential at N = 4 electrons per dot, because 
the electronic density distributions in the two cases are 
widely different. 

In order to include the e-e interactions in the study of 
AG structures we resort to density-functional theory [HI 
(DFT) within the 2D version of the local-density approxi- 
mation (2D-LDA) that has been shown to successfully de- 
scribe the electronic structure of individual and coupled 
QDs fabricated in the 2DEG [TH, [T3l - This takes us two 
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The electronic structure of the system was calculated 
with the Octopus package [22|. The code solves the 
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FIG. 1: (Color online) Left: Rectangular unit cell in real 
space. Right: Orthorhombic irreducible Brillouin zone (rect- 
angle FXSY) compared with the conventional one (YMK) 
for the honeycomb lattice. The basis vectors B x and B y are 
defined in the text. The hexagonal high symmetry points are 
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\~B y = K', M = T = M'. The mapping from the hexag- 



onal to the orthorombic cell is performed as follows: the Y — M 
line corresponds to the V — X — V path in the orthorhombic 
cell; the M — K segment corresponds to Y — K\ and the K — Y 
segment corresponds to the two segment path K' — K" — M\ 
with K" = | (B x + |B y ). Inset: part of the lattice in real 
space. 



steps beyond previous tight-binding studies. First, we are 
able to include the e-e interactions, thus producing the 
first self-consistent (density-dependent) band structure 
of AG. Second, we are also able to describe the system 
in a fully resolved manner in real space with a realistic 
model potential. Our results fully confirm the existence 
of Dirac points and massless Dirac fermions. However, 
we find that the threshold for the emergence of isolated 
Dirac points is increased in a way that significantly de- 
pends on the electron density. For example, from the 
non-interacting electron theory [tJ we know that both 
N = 1 and TV = 4 electrons per dot are candidates for 
the emergence of an isolated Dirac point, even though the 
N = 4 case requires a much deeper potential well. When 
e-e interactions are included we find that the threshold 
for the emergence of the Dirac point is slightly enhanced 
in the N = 1 case, but largely enhanced in the N = 4 
case. 

We consider electrons confined in GaAs/AlGaAs QDs 
in the effective mass approximation, which is the con- 
ventional approach to the modeling of individual and 
coupled semiconductor QDs [l6|, [l8j. We use the ef- 
fective mass m* = 0.067 mo and the dielectric constant 
e = 12.4 eo- Each QD is modeled by a cylindrical hard- 
wall potential with radius r = 52.5 nm and a tunable 
height Vb, and the lattice constant is fixed to a = 150 
nm (see Fig. [I]). These are values that can be reached 
experimentally flol-[2ll| and theyare similar to the previ- 
ous tight-binding study in Ref. [7] in order to allow direct 
comparison between the results. We point out that soft- 
ening the potential from a hard- wall to a Gaussian shape 
did not have qualitative effects on the results below. 



Kohn-Sham equations on a real-space discrete grid. The 
only required convergence parameter is the grid spacing, 
and it is therefore not bound to any particular choice of 
a basis set. In order to describe a truly 2D distribution 
of atoms we follow the procedure described in Ref. (23| . 
i.e., we impose a set of mixed boundary conditions (peri- 
odic in the plane, zero Dirichlet off-plane) and cut-off the 
Coulomb potential to zero along the direction perpendic- 
ular to the plane, while retaining its full long range of 
action within the plane. The levels obtained are rigor- 
ously equivalent to those calculated in an infinitely wide 
supercell in the direction perpendicular to the plane. 

As shown in Fig. [TJ the unit cell is chosen not to be 
minimal but to contain four dots, so that the 2D Bra- 
vais lattice becomes rectangular. The unit cell size is 
~ 130 nm x 225 nm and the grid spacing is ~ 2.45 nm. 
The reciprocal space cell is generated by the vectors 
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where a is the 
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interdot distance. The volume of this cell in reciprocal 
space is half the volume of the standard hexagonal BZ. 
Obviously, the physically meaningful results are not af- 
fected by our choice of the unit cell. However, in order to 
compare to calculations performed in the minimal (stan- 
dard) cell, our bands must be appropriately unfolded. 
The high symmetry lines of the standard BZ can all be 
mapped to corresponding paths into the smaller BZ. The 
mapping is described in the caption of Fig. [TJ For ease 
of comparison with previously published results, all the 
bands are displayed unfolded in this work. In the pro- 
cess of unfolding the band structure from the rectangular 
to the hexagonal cell, care must be exerted to avoid the 
phenomenon of aliasing, i.e., the spurious duplication of 
bands. This is greatly facilitated by the observation that 
the true bands must be continuous and differentiable at 
K along the Y — K — M line. 

As mentioned in the introduction, in order to assess 
the importance of e-e interactions, we compare the re- 
sults for noninteracting electrons with those computed 
using the Kohn-Sham DFT approach [15| within the 2D 
LDA. For the correlation part of the LDA we have used 
the parametrized form of the quantum Monte Carlo data 
calculated for the 2DEG by Attaccalite et al. [24| In view 
of previous works on QDs fabricated in the 2DEG HE|lZ| 
we believe that the 2D-LDA provides a reasonable ap- 
proximation for the energy bands considered here. 

In Fig. [2]we show the energy bands calculated for non- 
interacting (a) and interacting electrons in the LDA (b) 
when Vb = —0.6 meV. In both cases we find distinctive 
Dirac points at K with a linear dispersion relation. These 
are the defining attributes of graphene-like physics. From 
Fig. [2] it is clear that, in the case of singly occupied QDs 
in the AG, the e-e interactions do not make the Dirac 
point less stable. In general, the band structures are 
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FIG. 2: (Color online) Energy bands calculated for nonin- 
teracting (a) and interacting electrons (N = 1 per dot) in 
density-functional theory (b). Clear Dirac points and a linear 
dispersion relation can be found in both cases (circles). 
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FIG. 3: (Color online) Electron densities in artificial graphene 
(N — 1 per dot) calculated for two depths of the potential 
wells, V = -0.2 meV (up) and V = -0.5 meV (down). The 
left and right panels show the results without and with the 
e-e interactions, respectively. 



very similar, although, as expected, the e-e interactions 
reduce the bandwidth and thus increase the energy gaps 
between the two lowest bands. However, as an inter- 
esting exception to this tendency, the LDA result has a 
considerably smaller gap above the two lowest bands. We 
believe that this indicates the gap has a primarily kinetic 
origin, i.e. it arises from an overlap of neighboring local- 
ized orbitals forming a bonding-antibonding pair. Such 
a gap is reduced when the overlap of orbitals in the pair 
decreases, due to increased electron localization. 

The general tendency of increased electron localization 
due to e-e interactions is clearly visible in the electron 
density shown in Fig. [3l As expected, the relative dif- 
ference in the density with and without e-e interactions 




FIG. 4: (Color online) Energy gap at the M point in the 
Brillouin zone as a function of the quantum-dot depth. Re- 
sults without (squares) and with e-e interactions (circles) are 
shown for N = 1 per dot. The inset shows a zoom of the 
regime where the gap opens. 



is larger when \Vq\ is small. Indeed, for small \Vq\ the 
system is closer to the homogeneous 2DEG, where in- 
teractions (at small densities as in this case) drive the 
system to the Wigner crystal regime [25| . 

To obtain a closer view on the stability of the Dirac 
points and the effects of e-e interactions, we next ex- 
amine the onset and the size of the energy gap at the M 
point (see Fig. [2]). The gap is denned as the energy differ- 
ence between the Fermi level and the band closest to the 
Fermi level at the M point. Negative values correspond 
to crossings of the band(s) through the Fermi level, so 
that then there is no gap around the Dirac point. Fig- 
ure [31 shows the size of the gap as a function of |Vb|, i.e., 
the depth of the QDs. In the absence of interactions 
the threshold potential for the gap is Vo = —0.18 meV 
which agrees perfectly with the tight-binding result in 
Ref. [7]. When the interactions are included in the LDA 
level, the threshold shifts to Vq = —0.30 meV. In other 
words, the inclusion of e-e interactions leads us to pre- 
dict that a deeper QDs is required in order to achieve an 
isolated Dirac point in the system. When \Vq\ is further 
increased the results without and with interactions be- 
come very similar. The gap reaches its maximum value 
at Vo = —0.8 meV (or Vo = —0.6 meV in the absence 
of interactions), which can be regarded as the optimal 
QD depth for the stability of the Dirac point. As ex- 
pected, the gap shown in Fig. 0] closes asymptotically 
in the large- \Vo\ limit. The closing proceeds in a sim- 
ilar manner regardless of e-e interactions - albeit LDA 
results are not available at \Vo\ > 2 meV due to poor 
convergence. Gibertini et al argued that there is a lo- 
calization threshold in the system at \Vo\ ~ 3 meV, due 
to the formation of disorder-induced bound states. This 
effect is not included in our calculations, since the per- 
fect periodicity of the lattice keeps the wave functions 
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FIG. 5: (Color online) Energy bands close to the Fermi energy 
calculated for noninteracting (a) and interacting electrons, 
when there are four electrons in each quantum dot. The two 
lowest bands (not shown) are located at about —1.3 and —5.6 
meV in (a) and (b), respectively. 



extended, in principle up to the limit |Vb| — >• oo. 

Finally we consider occupying the QDs with more than 
one electron. As suggested in Ref. [7j, where a nonin- 
teracting tight-binding scheme was used, the next Dirac 
point should appear when TV = 4. In Fig. [5] we show 
the energy bands close to the Fermi energy when N = 4 
and Vo = —2 and -6.8 meV in the noninteracting (a) 
and interacting cases (b), respectively. The lowest two 
bands (not shown) are located at considerably lower en- 
ergies, i.e., at about —1.3 and —5.6 meV, respectively. 
It is noteworthy that although the Dirac point can be 
found in both cases, it is much less stable in the interact- 
ing system: first, due to the strong intradot interactions 
the QDs need to be significantly deeper than without 
interactions. Furthermore, the Dirac point in the inter- 
acting case appears in an energy gap that is smaller by 
about an order of magnitude. Consequently, decreasing 
the potential depth from \Vo\ = 6.8 meV to 6.5 meV al- 
ready disruptes the Dirac point, as the band above Ep 
in Fig. [5jb) is shifted downwards. In contrast, the Dirac 
point in the noninteracting case is stable down to about 
| Vo| = 1 meV. These results indicate that the realization 
of AG constitutes a challenge for the present engineering 
techniques if the QDs are occupied by several electrons. 

To summarize, we have studied the electronic prop- 
erties of artificial graphene fabricated in the two- 
dimensional electron gas. The electron-electron inter- 
actions, treated here within the two-dimensional local- 
density approximation in density-functional theory, gen- 
erally lead to stronger localization of the electrons in 
the quantum dots. The interactions shift the threshold 
for the emergence of isolated Dirac points to larger well 
depths than found without interactions. This effect is 
significantly pronounced when the number of electrons is 
increased. This sets a particular challenge for the realiza- 
tion of Dirac points if the quantum dots in the hexagonal 
lattice are occupied by more than a single electron. 

This work has been supported by the Academy of Fin- 
land, the Wihuri Foundation, and the DOE grant DE- 
FG02-05ER46203 (SP). CSC Scientific Computing Ltd. 
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